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The term “surface polarization” is introduced to describe the in-plane polarization existing at the 
surface of an insulating crystal when the in-plane surface inversion symmetry is broken. Here, the 
surface polarization is formulated in terms of a Berry phase, with the hybrid Wannier representation 
providing a natural basis for study of this effect. Tight binding models are used to demonstrate 
how the surface polarization reveals itself via the accumulation of charges at the corners/edges for 
a two dimensional rectangular lattice and for GaAs. 
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For over two decades, it has been understood that 
the electric polarization P of an insulating crystal is 
a bulk quantity whose electronic contribution is de¬ 
termined modulo 2eR/^2 (where R is a lattice vector 
and Q is the unit cell volume) by the Bloch functions 
through a Berry-phase expression, or alternatively, in real 
space through the charge centers of the Wannier func¬ 
tions BE]. It was also shown that the macroscopic sur¬ 
face charge of an insulating crystal is predicted by the 
standard bound-charge expression = P • n (where n 
is the surface normal) |3], as illustrated schematically in 
Fig. Ha). 

Here, we introduce and analyze a related quantity, the 
‘^surface polarization, ” defined as a 2-vector P lying in 
the plane of an insulating surface of an insulating crystal. 
By analogy with the bulk 3-vector P, it has the property 
that when two facets meet, the linear bound-charge den¬ 
sity appearing on the shared edge is predicted to be 

A^dge = Pi . ni + P2 ■ n2 (1) 

where Pj is the surface polarization on facet j and is 
a unit vector lying in the plane of the facet and point¬ 
ing toward (and normal to) the edge, as illustrated in 

Fig-Hb). 

This surface polarization P is quite distinct from the 
dipole per unit area normal to the surface, which has also 
been called “surface polarization” by other authors m- 
The latter is always present regardless of the symmetry of 
the surface, and manifests itself macroscopically through 
the surface work function. In contrast, our surface po¬ 
larization P lies in-plane and is nonzero only when the 
symmetry of the terminating surface supports a nonzero 
in-plane vector, as for example on the (110) surface of 
GaAs. It can also arise from a spontaneous symmetry¬ 
lowering surface reconstruction, as observed recently at 
the Pbi_a,Sna,Se (110) surface [6] and predicted for an 
ultrathin film of SrCrOs on SrTiOs substrate (001) [7]. 
The surface polarization will be most evident when the 
bulk P vanishes, as will be the case for the systems dis¬ 
cussed below. 

The purpose of this Letter is to extend the Berry-phase 
theory to the case of surface polarization P as defined 



FIG. 1. Illustration of bulk and surface polarization effects. 
The polarizations are denoted by black arrows, and net pos¬ 
itive and negative bounded charged are in red and blue, re¬ 
spectively. (a) Bulk polarization gives rise to surface charges 
a. (b) Surface polarization gives rise to edge charges A. 


above. To do this, we introduce a formulation based on 
hybrid Wannier functions (HWFs), which are Bloch-like 
parallel to the surface and Wannier-like in the surface- 
normal direction IHHII]- This allows for the use of Berry- 
phase methods parallel to the surface while allowing a 
real-space identification of the surface-specific contribu¬ 
tion in the normal direction. We illustrate the concept 
first for a “toy” 2D tightbinding (TB) model, demon¬ 
strating the method of calculating the surface polariza¬ 
tion. We then consider an atomistic 3D model of an ideal 
(110) surface of a generic III-V zincblende semiconduc¬ 
tor, using a TB model of GaAs to describe the electronic 
structure. In both cases, we confirm that the surface 
polarization correctly predicts corner and edge charges. 

We first show how to express the surface polarization 
in terms of the Berry phases of HWFs for a 2D insulat¬ 
ing sample, which we take to lie in the (x, z) plane. We 
take the “surface” (here really an edge) to be normal to i 
and introduce HWFs \hij{kx)), where I indexes unit-cell 
layers normal to the z direction and j runs over occupied 
Wannier functions in a single unit cell. For the bulk, 
the lattice is periodic in 2 : as well as x, and the \hij{kx)) 
and their centers z]^^^^{kx) can be obtained using the ID 
construction procedure given in Ref. [12]. To study the 
surface behavior we consider a ribbon consisting of a fi¬ 
nite number of unit cells along z. We then construct and 
diagonalize the matrix Z^n = {i^m{kx)\z\2pri{kx)), whose 
eigenvectors yield the HWFs and whose eigenvalues give 
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the HWF centers zij{kx). In practice these are easily 
identified with the bulk covering the range of I 

values that define the ribbon, with only modest shifts in¬ 
duced by the presence of the surface, allowing a common 
labeling scheme for both. 

If we were interested in computing the dipole moment 
normal to the surface, we could obtain this from an anal¬ 
ysis of the /Ca^-averaged ^ positions zij of the HWFs, where 
a is the lattice constant along x. However, our purpose 
here is different: we want to compute the polarization 
parallel to the surface. For this, we compute the Berry 
phase 


'yx,ij dkx l^b)’ (^) 

of each HWF “band” (Ij) as kx runs across the ID BZ. 
Doing the same for the bulk HWFs (these are indepen¬ 
dent of 1) and taking the difference, we obtain a set of 
Berry-phase shifts Ajx,ij = lx,ij — from which the 

electronic surface polarization can be determined via 

top surf 

= E (3) 

Zjf=center 

where a factor of two has been included for spin degener¬ 
acy and O is the edge repeat length a in 2D. Since A^x,ij 
decays exponentially into the bulk, the sum will converge 
within a few layers of the surface, but for definiteness we 
sum to the center of the ribbon. If the zij values of some 
neighboring HWF bands overlap, the procedure needs to 
be generalized by grouping the HWFs into layers and us¬ 
ing a multiband generalization to assign contributions to 
each layer. 

The generalization to a 3D crystal with surface nor¬ 
mal to 2 : is straightforward. The HWFs are \hij{kx^ky)) 
with centers zij{kx^ ky). The surface polarization is 
then obtained by computing Berry phases with respect 
to kx as before, averaging over all ky^ and multiplying by 
the lattice constant a divided by the surface cell area O. 
The other surface polarization Vy is given by the same 
formalism but with x and y reversed. 

In the models considered in this paper, the surface 
polarization is purely electronic, as the ions are held 
fixed in their bulk positions. More generally, Vx = 
with the ionic contribution given by = 
Eir Zr (Xir - where and Xi^ are the x 

position and bare charge of ion r in cell /, and is 

the corresponding bulk position of the same atom. 

To illustrate these ideas, we start by considering a 
tight-binding (TB) model of the simple 2D crystal shown 
in Fig.j^a). We assume a rectangular lattice with an as¬ 
pect ratio 6/a = 0.8. There are two atoms symmetrically 
located along a diagonal of the rectangular unit cell with 
coordinates (|, |) and (|,|), so that the bulk crystal has 
inversion symmetry. We consider only one s orbital per 
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(d) 



FIG. 2. (a) Illustration of the TB model, where four unit cells 
are presented. The atoms are denoted by black dots. Nearest 
neighbor hoppings ti, t 2 and ts are shown in solid blue lines. 
Next-nearest neighbor hoppings and ts are shown in dashed 
blue lines. V are shown by black arrows. The induced Qcomer 
are denoted by red (positive) and blue (negative) large dots 
at the corners, (b) Band structure of the TB model in the 
{kx,kz) space. The inset shows the high symmetry points in 
the 2D Brillouin zone, where F, X, M, X' refer to (0,0), (0,|), 
(|,|) and (|,0), respectively, (c) Band structure along kx for 
the 2D slab model that is infinite along x while 10-cell-thick 
in z. (d) Difference between effective x positions of each HWF 
and that deep in the bulk. 


atom with onsite energy taken to be zero, and assume 
that each atom contributes one electron so that only the 
lower band is (doubly) occupied. We take the nonzero 
hoppings to be those shown in Fig.[^a) and choose their 
values to be ti = —2.2, ^2 = —0.15, and = —0.1, 
^4 = —0.09 and = —0.06 in eV. The position operators 
are taken to be diagonal in the local-orbital representa¬ 
tion so that {(l)i\z\(l)j) = Zidij. 

We plot the bulk band structure of the TB Hamilto¬ 
nian in Fig. [^b). For the selected parameters the band 
gap is large compared to the band widths; in particular, 
the upper (unoccupied) band is quite fiat. Next we com¬ 
pute the surface polarization of a ribbon cut from the 
2D lattice, taking it to be ten unit cells thick along z 
and infinite along x. For the atoms in the surface lay¬ 
ers, the hoppings to the interior atoms are the same as 
those described above, while the hoppings to the vac¬ 
uum side are set to zero. We used an equally spaced 
60-point kx grid. At each kx the 20 x 20 Hamiltonian 
is diagonalized, resulting in the band structure shown 
in Fig. [^c). There are no obvious surface states, and 
in fact the result looks almost indistinguishable from a 
surface projection of the bulk band structure. The eigen¬ 
functions I'lpnikx)) are expressed in the tight-binding ba- 
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sis as \ipn{k^)) = Y.jCnj{kx)\Xj{kx)), where the \Xj{kx)) 
are the Bloch basis function formed as a Fourier sum 
at wavevector kx of atomic orbitals |0i). From the ten 
occupied bands we construct the 10 x 10 position ma- 
trix Zmn = {i’m{kx)\z\'lpn{kx)) = ZjC*^j{kx)Cnj{kx) ■ 
Diagonalizing this matrix, we get ten eigenvalues zi{kx) 
that can each be clearly associated with a particular unit 
cell layer, and ten eigenfunctions that are the HWFs. 
We label the HWF \hi{kx)), where I is the layer index 
running from the bottom to the top of the ribbon. 

Next we calculate jx,h fhe Berry phase along kx^ for 
each I using Eq. §. Deep in the interior these Berry 
phases become equal to tt within numerical precision, 
while the Berry phases near the edge are slightly shifted 
away from tt, leading to a nonzero surface polarization as 
shown in Fig. [^d). 

The value of the surface polarization obtained from 
Eq. ([^ is = ±2.1 x 10“^ e for the top and bottom 
surfaces respectively. Similarly we can compute the sur¬ 
face polarizations for the left/right surfaces using a rib¬ 
bon ten cells wide in x and infinite along z. We obtain 
Vz = ±4.7 X 10“^ e along the left and right edges respec¬ 
tively. At the corners, the surface polarizations are di¬ 
rected head-to-head or tail-to-tail, as shown in Eig.j^a). 

Given the values of the surface polarizations in the 2D 
model, we predict that the charge accumulation at the 
corner of a finite sample should equal the sum of the 
two adjacent surface polarizations, here \Vx\ I'Pz] = 
6.8 X 10“^ e. To test this, we directly calculate the cor¬ 
ner charge in a finite 2D sample, specifically a 10 x 10 
supercell, large enough to ensure neutrality in the cen¬ 
tral region and in the middle of the edges of the sam¬ 
ple. The corner charge is obtained by summing up 
the on-site charge differences, relative to the bulk, for 
atoms in the quadrant containing the corner. We find 
Qcorner = 6.8 X 10“^ 6 for the top left and bottom right 
corners, and —6.8 x 10“^ e for the other two corners, in 
agreement with our prediction from the computed surface 
polarizations. 

We now consider a TB model of a generic III-V 
zincblende semiconductor, with GaAs as the prototypical 
example. The crystal structure is characterized by Ga-As 
zigzag chains running along (110). Although the crystal 
structure does not have inversion symmetry, the tetra¬ 
hedral symmetry forbids a nonzero spontaneous polar¬ 
ization. We use tight-binding parameters from Ref. [13], 
in which is shown the bulk bandstructure and density of 
states. The unit cell contains two Ga and two As atoms, 
each with four sp^ hybridized orbitals and four electrons, 
as shown in Eig. [^a). The position matrix is assumed 
to be diagonal and atom-centered in the basis of tight- 
binding orbitals. la 

To describe the (110) surface, we construct a slab ge¬ 
ometry as shown in Eig. [^a), and we henceforth label 
the Gartesian directions as shown there. That is, the 
surface, which is normal to i, has zigzag chains running 



FIG. 3. (a) Illustration of the GaAs slab studied in the TB 
model, where the blue and green balls represent Ga and As 
atoms, respectively. The grey shaded planes denote the (110) 
family planes, (b) Electronic band structure of the GaAs slab 
in the 2D Brillouin zone, with the thickness of 8 cells z. (c) 
Difference between the 'jx of each group of HWFs and that 
deep in the bulk. 


along y. Since the two atoms making up these chains are 
inequivalent, we expect a surface polarization in the x 
direction. We take the slab to be eight unit cells thick; 
for the atoms in the surface layer, the hoppings to the 
atoms inside the slab are the same as in the bulk, while 
the hoppings to the vacuum side are set to zero. At each 
{kx^ky) of the 100 x 100 k grid in the surface BZ, the 
128 X 128 Hamiltonian is diagonalized, and we obtain the 
band structure for the slab, shown in Fig. [^b). Surface 
states are evident as isolated bands. 

Next, we diagonalize the 64x64 position matrices 
Z{kx^ky) constructed from the eigenstates of the occu¬ 
pied bands. The eigenvalues, which are the ^ coordinates 
of the HWF centers, can be clearly divided into groups, 
each consisting of four HWFs j representing the four Ga¬ 
As bonds around an As atom, each group being associ¬ 
ated with one of the 16 atomic layers 1. In this case, it 
is more useful to calculate the Berry phase of each group 
of HWFs rather than of each single HWF |3| . 

As expected, the Berry phase in the y direction along 
the zigzag chain is found to be zero, but in the x direction 
it is nonzero for the HWF groups near the top and bot¬ 
tom surfaces of the slab. Thus, we confirm that there is 
a nonzero surface polarization Vx- We plot the difference 
between the Berry phase ^x of each group of HWFs and 
that for the bulk in Fig.j^c). By summing up the contri- 
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FIG. 4. (a) Overhead view of the hexagonal GaAs nanowire. 
The dashed black line, which meets the edges along [111] at 
an angle of ^ = 35.26°, shows the direction along which the 
zigzag surface chains run. The relevant surface polarizations 
at the side surfaces are denoted by black arrows. The blue 
and red vertical edges mean net negative and positive edge 
charge distributions, respectively, (b) On-site charge distri¬ 
bution summed over the trilayer. Red and blue dots represent 
positive and negative net charges, respectively. The sizes of 
dots indicate the magnitudes of the on-site net charge, the 
left and right regions to the dashed vertical line show the total 
and symmetric part of onsite charges, respectively. 


butions from each group of HWFs from the center of the 
bulk to one surface, the total surface polarization is found 
to be 0.178 ejL. Here L = a/\/2 is the repeat length of 
the zigzag chain, i.e., the surface cell dimension along 
where a is the surface lattice constant along x. Subdi¬ 
viding the dominant surface-group contribution further, 
we find that the surface polarization comes mainly from 
the surface-most HWF, corresponding to a shift of the 
center of the dangling bond on the surface As atom. 

The surface polarization on the {110} surfaces predicts 
an accumulated line charge for the common edge of two 
such surfaces. In order to demonstrate this effect, we con¬ 
sider a hexagonal wire of GaAs that is infinite along [111], 
with a periodicity corresponding to three of the GaAs 
buckled (111) layers. In this case, the six side surfaces of 
the wire are all {110} planes: (110), (101), (Oil), (110), 
(lOl), and (oil). As shown in Fig.a), on each side facet 
the surface polarization is perpendicular to the zigzag 
chains, forming a pattern of V vectors shown as black ar¬ 
rows. The surface polarizations for each neighboring pair 
of side facets have a common component along [111], but 
are head-to-head or tail-to-tail for the component normal 
to [111], leading to alternating positive and negative line 
charges for the six edges as shown. According to Eq. 
we expect the line charge per three-layer vertical period 
to be Qsl = 27^ cos 6>-3L cos 0 = 0.71e, where the 3L cos 0 
factor is the vertical period. 

For comparison, we directly calculate the edge charges 
per trilayer period in a nanowire with a radius of 8 atoms. 


We sum up the site populations within the TB model 
with a 60-point k grid along [111]. The onsite charge is 
the difference from the bulk value. The computed onsite 
charges are shown in the left half of Fig. |^b), while the 
right half shows the corresponding results after averaging 
with a 60°-rotated version of itself. The surface charges 
decay rapidly into the bulk, leading to a neutral bulk 
state inside the nanowire. Also, a surface dipole density 
normal to the surface is clearly visible, especially in the 
orientationally averaged results. However, we are inter¬ 
ested in the accumulation at the edges, which is obviously 
present in the unaveraged results in the left half of the 
figure. The edge charge is calculated by summing up the 
onsite charges in the wedge-shaped region illustrated in 
Fig. gb), using a weight of 1/2 for atoms located on its 
radial edges. The total edge charge per trilayer is found 
to be ±0.71e, in agreement with the value predicted using 
the previously calculated surface polarization. 

We emphasize that this numerical value is not intended 
to be realistic for GaAs. A more accurate estimate would 
require the use of an improved tight-binding model and 
treatment of surface relaxations and dielectric screening 
effects, or better, direct first-principles calculations. Our 
purpose here has been to show that the surface polariza¬ 
tion as defined here correctly predicts edge charges. We 
note that an analysis based on maximally localized Wan- 
nier functions na is also possible. However, we believe 
our HWF-based approach is more natural, as the Wan- 
nier transformation is only done in the needed direction 
and no iterative construction is required. 

We stress that the concept of surface polarization V 
is quite general, occurring whenever the surface symme¬ 
try is low enough. In some cases this can arise from 
a spontaneous symmetry-lowering surface relaxation or 
reconstruction, allowing “surface ferroelectricity” if it is 
switchable. In other cases, as for GaAs (110), the ideal 
surface space group already has low enough symmetry 
to allow a nonzero P. This will occur quite generally 
for low-angle vicinal surfaces. The concept also applies 
to planar defects such as domain walls, stacking faults, 
and twin boundaries, and to heterointerfaces; if V is 
present within this plane, it may induce a line charge 
where the plane intersects the surface. Such edge and 
line charges are potentially observable using electric force 
microscopy m , electron holography m, or other ex¬ 
perimental methods. Finally we note that the concept 
of surface polarization may become more subtle in the 
presence of orbital magnetization, which we have omit¬ 
ted from our considerations here. 

In summary, we have formulated the concept of surface 
polarization, i.e., the dipole moment per unit area paral¬ 
lel to the surface, which can exist whenever the surface 
symmetry is low enough. Using TB models we have com¬ 
puted the surface polarizations for a 2D toy model and a 
generic HI-V zincblende semiconductor, and shown that 
the predicted corner or edge charges are in good agree- 
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ment with direct calculations. We point out that surface 
and interface polarizations can be responsible for observ¬ 
able effects, and perhaps even desirable functionality, in 
a broad range of insulating materials systems. 
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